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ABSTRACT 

Spectroscopic observations of the low-mass eclipsing binary NSVS 11868841 
have been obtained and the radial velocities were derived for both components. 
The masses and radii determined for the components are Mi =O.87O±O.O74M0, 
M2=0.607±0.053MQand Ri =0.983 ±0.030 Rq, i?2 =0.901±0.026 R©. Both 
the primary and secondary stars' radii are 10 % and 57% larger than those 
i-^ ' of zero-age-main-sequence stars with the same masses. This discrepancy may 

be arisen from the large spot coverage of both stars. We collected abso- 
O I lute parameters of 21 low mass double-lined eclipsing binaries and compared 

-*— * ■ their positions in the mass-radius and mass-effective temperature panels. The 

C^ , large radii and lower effective temperatures are solved neither with difference 

in metallicity nor in mixing length parameters. These discrepancies in the low 

^SJ ■ mass stars may be originated by magnetic fields causing inhibition convective 

^ ' energy transport which leads to large magnetic spot coverage on the surface 

^. of a low mass star. 

^^ ■ 

\f\ • Key Words: Stars; activity — Stars: fundamental parameters — Stars: low- 

I ^ \ mass — Stars: eclipsing 

o. 

O . 1. INTRODUCTION 



Eclipsing binaries are the best laboratories for determining the fundamen- 
tal physical properties of stars. Detached, double-lined eclipsing binaries yield 
\^ ' direct and accurate measures of the masses and radii of the component stars. 



X 



C^ 



Measuring of these quantities has always been observationally challenging. 
Large-scale surveys, providing imaging and photometric catalogues for sub- 
stantial fractions of the celestial sphere, are now available at a wide range 
of wavelengths. Over the last few years, these surveys have become invalu- 
able research tools for investigating the properties of intrinsically rare objects. 
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These surveys have proven most useful in searches for low-mass stars. Stars 
with masses less than 1 M0 are regarded as low-mass stars and have late spec- 
tral types, namely G, K or M. We cannot vary the conditions prevalent on an 
individual star, so observations of many individual stars are combined to build 
up a picture of how stars of different mass and composition evolve. Since the 
intermediate- and high-mass stars are intrinsically bright sources most studies 
have centered on them with little consideration of M dwarfs with masses below 
0.6 Mq. This reflects the availability of more precise observational constraints, 
and also greater analytic tractability in modelling higher-mass stars. Recently, 
however, the lower main-sequence has attracted more attention, with a series 
of detailed models extending past the hydrogen-burning limit to the bound- 
ary between low-mass brown dwarfs and giant planets. Understanding the 
structure and evolution of stars is a basic goal of stellar astronomy. Critical 
tests of the evolution theory for stars can be made on a small set of eclips- 
ing binary stars (Baraffe et al. 1998, Girardi et al. 2002, Siess et al. 1997, 
Chabrier et al. 2007). While the results of stellar evolution models compare 
favourably to data for main-sequence stars with masses greater than that of 
the Sun, evidence has been growing that the models for stars on the lower 
main-sequence have problems when confronted with precise masses and radii 
from double-lined eclipsing binaries. As shown by Torres and Ribas (2002) 
the evolutionary models underestimate the radii and overestimate the effec- 
tive temperatures of the low-mass stars. These discrepancies were confirmed 
in subsequent studies by Ribas(2003), Lopez-Morales and Ribas (2005), Tor- 
res et al. (2006), Morales et al. (2009a), Morales et al. (2009b), and others. 
Hints of the radius discrepancies have also been described by Popper (1997), 
Clausen et al. (1999), and others. Recently, MuUan and MacDonald (2001), 
Torres et al. (2006), Morales, Ribas and Jordi (2008) proposed a hypothesis 
based upon the effects of stellar activity to explain the discrepancies. Fur- 
thermore Stassun et al. (2009) called attention to that low-mass stars appear 
older or younger in mass-radius diagram depending on whether post- or pre- 
main-sequence models are used. As the number of low-mass eclipsing binaries 
studied carefully, both photometrically and spectroscopically, is increased a 
better comparison with the theoretical models can be made. 

Light variability of NSVS 11868841 ( a = 23''17™58% d = + 19°17'03"(J2000), 
hereafter NSVS 1186) was detected by Wozniak et al.(2004) from the North- 
ern Sky Variability Survey. They classified it as a detached eclipsing 
binary with combined median magnitude of V=14™.084±0.076 and an or- 
bital period of 0.60179 days. High-precision light curves in the Johnson V-, 
R-, and I-bandpass were obtained by Coughlin and Shaw (2007). Since the 
radial velocities of the components were not available they have determined 
the radii of the stars by analysing the V- and R-bandpass light curves. In 
this study we present the first radial velocities for the components of double- 
lined, short-period eclipsing binary NSVS 1186 and the physical parameters 
obtained from the analyses of radial velocities and light curves. Both com- 
ponents show high level chromospheric activity, evidenced by H^ emission 
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lines in their spectra. We discuss the properties of the binary and in rela- 
tion to the discrepancies between empirical and theoretical mass-radius and 
mass-effective temperature relations for low-mass stars. 

2. SPECTROSCOPIC OBSERVATIONS 
2.1. Spectroscopy 

Since the wide-band light curves of NSVS 1186 are available we observed 
the star spectroscopically. Optical spectroscopic observations were obtained 
with the Turkish Faint Object Spectrograph Camera (TFOSC) attached to 
the 1.5 m telescope on August 23, 24, and 25, 2008 under good seeing con- 
ditions. Further details on the telescope and the spectrograph can be found 
at http : //www .tug . tubitak . gov . tr'. The wavelength coverage of each spec- 
trum was 4100-8100 A in 11 orders, with a resolving power of A/AA 7000 
at 6563 A and an average signal-to-noise ratio (S/N) was ~150. We also 
obtained a high S/N spectrum of the M dwarf GJ 740 (MO V) and G J 623 
(Ml. 5 V) for use as templates in derivation of the radial velocities (Nidever 
et al. 2002). 

The electronic bias was removed from each image and we used the 'crre- 
ject' option for cosmic ray removal. The echelle spectra were extracted and 
wavelength calibrated by using a Fe-Ar lamp source with help of the IRAF 
ECHELLE package. The stability of the instrument was checked by cross cor- 
relating the spectra of the standard star against each other using the FXCOR 
task in IRAF. The standard deviation of the differences between the velocities 
measured using FXCOR and the velocities in Nidever et al. (2002) was about 
1.1 km s^^. 

2.1.1. Spectral Types and Temperature Estimates 

We have used the spectra to reveal the spectral type of the primary compo- 
nent of NSVS 1186. For this purpose we have degraded the spectral resolution 
from 7 000 to 3 000, in order to use the calibrations by Hernandez et al. (2004) 
by convolving them with a Gaussian kernel of the appropriate width, and we 
have measured the equivalent widths (EW) of photospheric absorption lines 
for the spectral classification. We have followed the procedures described by 
Hernandez et al., choosing metallic lines in the blue- wavelength region, where 
the contribution of the primary component is considerably larger than that 
of the secondary in the spectra. From several spectra we measured equivalent 
widths of some spectral lines assuming the light contribution of the primary 
star is 0.6, as given in Table 1. 

From the calibration relations of i^W^-Spectral-type given by Hernandez 
et al. (2004), we have derived a spectral type of G9±l for the primary com- 
ponent. The effective temperature deduced from the calibrations of Drilling 
& Landoh (2000), de Jager & Nieuwenhuijzen (1987) and Popper (1980) is 
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TABLE 1 

EQUIVALENT WIDTHS OF THE SELECTED LINES IN THE 

SPECTRA. 

Spectral lines EWp (A) EW^ (A) 

Ca I A 4226 2.13±0.11 1.39±0.12 

Fe I A 5329 l.lliO.ll 0.99±0.12 

Fe I A 4271 1.44±0.13 1.41±0.12 

Ca i+Fe I A 5270 2.91±0.21 2.11±0.21 

Ca I A 6162 1.99±0.42 1.11±0.09 

Fe I A 5079 0.79±0.07 0.47±0.11 



5 230±80, 5 230±80, and 5 330±141K for the primary component, respec- 
tively. The mean effective temperature of the primary component deduced 
from the spectra is 5 250±95 K . 

NSVS1186 is Usted in several large photometric databases consolidated 
in the Naval Observatory Merged Astronomical Dataset (NOMAD- 
1.0, Zacharias et al. 2004), which provide optical magnitudes of B=14™.77±0™.02, 
V=13'".87±0'".01, R=13"'.94±0'".G1, 1=12™. 74±0™.02. Since the magnitudes 
collected from photographic measurements and the colors are inconsistent no 
attempt has been made to calculate the effective temperatures of the compo- 
nents. The infrared magnitudes are taken from 2MASS (Cutri et al. 2003) cat- 
alog as J=12'".518±0™.023, H=12™.069±0™.021, and K=ll'".995±0™.019. 
The observed infrared colours of J-H=0™.449±0'".031 and H-K=0".074±0™.028 
correspond to a combined spectral type of G9±2 is in a good agreement with 
that we derived by spectral lines alone. Using the depth of the eclipses in three 
bands we estimate light contribution of the primary component as 0.60, 0.57, 
0.55 and 0.56 in the V, J, H and K bands respectively. Hence V-K, J-H and 
H-K colors of the primary component corresponds to a spectral type of G9±2. 
We estimated a temperature of 5 240 ±230 K from the calibrations of Toku- 
naga (2000). Temperature uncertainty of the primary component results from 
considerations of spectral type uncertainties, and calibration differences. The 
effective temperature of the primary star what we derived from the infrared 
photometric measurements is in a good agreement with that we estimated 
from the spectra alone. The weighted mean of the effective temperature of 
the primary component is 5 250±135K. 

2.1.2. Rotational Velocities of the Components 

The width of the cross-correlation profile is a good tool for the measure- 
ment of usini (see, e.g., Queloz et al. 1998). The rotational velocities {vsmi) 
of the two components were obtained by measuring the FWHM of the CCF 
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peaks in nine high-S/N spectra of NSVS 1186 acquired close to the quadra- 
tures, where the spectral lines have the largest Doppler-shifts. In order to 
construct a cahbration curve FWHM-i;sinj, we have used an average spec- 
trum of HD 27962, acquired with the same instrumentation. Since the rota- 
tional velocity of HD 27962 is very low but not zero (usini ~11 km s~^, e.g., 
Royer, Zorec & Fremat 2004 and references therein), it could be considered 
as a useful template for low-mass stars rotating faster than vsini ~ 10 km 
s~^. The spectrum of HD 27962 was synthetically broadened by convolution 
with rotational profiles of increasing v sin i in steps of 5 km s^^ and the cross- 
correlation with the original one was performed at each step. The FWHM 
of the CCF peak was measured and the FWHM-w sini calibration was estab- 
lished. The usini values of the two components of NSVS 1186 were derived 
from the FWHM of their CCF peaks and the aforementioned calibration re- 
lations, for few wavelength regions and for the best spectra. As a result the 
projected rotational velocity of 77±3 km s^^ for the primary star and 62±4 
km s~^ for the secondary star were estimated. However, note that the spec- 
tral lines of the components could not be separated even at the quadratures. 
Thus the projected rotational velocities estimated in this way should involve 
higher uncertainties than those given here. 

3. ANALYSIS 

3.1. The Orbital Configuration 

To derive the radial velocities for the components of binary system, the 
16 TFOSC spectra of the eclipsing binary were cross-correlated against the 
spectrum of GJ 740, a single-lined MOV star, on an order-by-order basis us- 
ing the FXCOR package in IRAIO- The majority of the spectra showed two 
distinct cross-correlation peaks in the quadrature, one for each component of 
the binary. Thus, both peaks were fit independently in the quadrature with a 
Gaussian profile to measure the velocity and errors of the individual compo- 
nents. If the two peaks appear blended, a double Gaussian was applied to the 
combined profile using de-blend function in the task. For each of the 16 obser- 
vations we then determined a weighted-average radial velocity for each star 
from all orders without significant contamination by telluric absorption fea- 
tures. Here we used as weights the inverse of the variance of the radial velocity 
measurements in each order, as reported by FXCOR. In these data, we find no 
evidence for a third component, since the cross-correlation function showed 
only two distinct peaks. We adopted a two-Gaussian fit algorithm to resolve 
cross-correlation peaks near the first and second quadratures when spectral 
lines are visible separately. Figure 1 shows examples of cross-correlations ob- 
tained by using the largest FWHM at nearly first and second quadratures. 
The two peaks, non-blended, correspond to each component of NSVS 1186. 



^IRAF is distributed by the National Optical Observatory, which is operated by the 
Association of the Universities for Research in Astronomy, inc. (AURA) under cooperative 
agreement with the National Science Foundation 
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Fig. 1. Sample of Cross Correlation Functions (CCFs) between NSVS1186 and the 
radial velocity template spectrum around the first and second quadrature. 



The stronger peaks in each CCF correspond to the more luminous component 
which has a larger weight into the observed spectrum. 

The heliocentric radial velocities for the primary (Vp) and the secondary 
(Vs) components arc listed in Table 2, along with the dates of observation 
and the corresponding orbital phases computed with the ephemeris given by 
Coughlin and Shaw (2007). The velocities in this table have been corrected 
to the heliocentric reference system by adopting a radial velocity of 9.5 km 
s~^ for the template star GJ740. The radial velocities listed in Table 2 are 
the weighted averages of the values obtained from the cross-correlation of 
orders #4, #5, #6 and #7 of the target spectra with the corresponding order 
of the standard star's spectrum. The weight Wi = l/cr^ has been given to 
each measurement. The standard errors of the weighted means have been 
calculated on the basis of the errors (ui) in the radial velocity values for each 
order according to the usual formula (e.g. Topping 1972). The Ui values are 
computed by FXCOR according to the fitted peak height, as described by Tonry 
& Davis (1979). 

First we analysed the radial velocities given in Table 2 for the initial or- 
bital parameters. We used the orbital period held fixed and computed the 



THE ABSOLUTE PARAMETERS FOR NSVS 11868841 7 

TABLE 2 

HELIOCENTRIC RADIAL VELOCITIES OF NSVS 1186. THE 

COLUMNS GIVE THE HELIOCENTRIC JULIAN DATE, THE ORBITAL 

PHASE, THE RADIAL VELOCITIES OF THE TWO COMPONENTS 

WITH THE CORRESPONDING STANDARD DEVIATIONS. 



HJD 2400000+ 


Phase 


Star 1 




Star 2 








V, 


a 


Vs 


a 


55069.4412* 


0.0536 


-120.0 


14.3 






55069.4855 


0.1272 


-162.0 


9.2 






55069.5284 


0.1985 


-193.6 


4.3 


80.6 


4.2 


55069.5727 


0.2721 


-198.1 


3.9 


85.6 


6.3 


55068.3810 


0.2918 


-183.7 


3.2 


80.0 


8.4 


55068.4258 


0.3663 


-171.6 


3.0 


47.2 


8.6 


55068.4707 


0.4409 


-128.7 


13.2 






55068.5134 


0.5119 


-68.2 


14.9 






55068.5578 


0.5857 


-28.6 


11.1 






55067.3674 


0.6076 


-6.0 


7.3 


-184.6 


19.7 


55068.6006 


0.6568 


19.2 


3.3 


-214.3 


12.6 


55067.4122 


0.6820 


28.3 


3.1 


-235.1 


10.2 


55067.4552 


0.7534 


40.0 


1.4 


-244.3 


13.3 


55067.4976 


0.8240 


19.9 


9.2 


-230.6 


5.2 


55067.5419 


0.8975 


1.1 


2.4 






55067.5846* 


0.9684 


-56.5 


15.5 







* These data were not used in the RV analysis. 



eccentricity of the orbit, systemic velocity and semi-amplitudes of the radial 
velocities. The results of the analysis performed with PHOEBE are as follows: 
e=0.002±0.001, i.e. formally consistent with a circular orbit, 7=-80±3 km 
s^^, iiri = 118±5 and iir2=169±6 km s~^. Using these values we estimate the 
projected orbital semi-major axis and mass ratio as: asini =3.412±0.093 
R© and q = ^=0.70±0.04. The observed and calculated radial velocities of 
both components are plotted in Figure 2 for comparison. 



3.2. Modelling the Light Curves 

As we mentioned in § 1 the V-, R- and I-bandpass light curves were ob- 
tained by Coughlin and Shaw (2007). The differential magnitudes and orbital 
phases were read out carefully by the method of dexteeo. We used the most 
recent version of the eclipsing binary light curve modelling algorithm of Wil- 
son & Devinney (1971, with updates), as implemented in the phoebe code 



^http : / /Ac . zah ■ unl-heidelberg . de/sdext er ' 
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Fig. 2. Radial velocity curve folded on an orbital period of 0.601790 days, where 
phase zero is defined to be at primary mid-eclipse. Symbols with error bars show the 
RV measurements for two components of the system (primary: circles, secondary: 
squares). Open circles were excluded from the RV analysis. 



of Prsa & Zwitter (2005). The code needs some input parameters, which de- 
pend upon the physical structures of the component stars. In the light curve 
solution we fixed some parameters whose values can be estimated from global 
stellar properties, such as effective temperature of the primary component 
and mass-ratio. The weighted mean of the effective temperature has been 
estimated to be 5 250±135K and the mass-ratio of 0.70±0.04. Using these 
values we adopted the linear limb-darkening coefficients from Van Hamme 
(1993) as 0.39 and 0.28 for the primary and secondary components, respec- 
tively; the bolometric albedos from Lucy (1967) as 0.5, typical for a convective 
stellar envelope, the gravity brightening coefficients as 0.32 for the both com- 
ponents. The rotation of components is assumed to be synchronous with the 
orbital one. We start with the adjustable parameters: the orbital inclination, 
i, the effective temperature of the less massive star, Tg/^^j the potential of 
the components, ili and Q2, and the monochromatic luminosity of the more 
massive star, Li. Assuming it is a detached binary we used Mode — 2 of the 
DC program for individual solution. First we determined the preliminary ele- 
ments and compared with the observed light curves. Since the observed light 
curves show distortions we started with the two spot solution, one on the 
primary and other on the secondary component. The initial spot parameters 
were adopted from Coughlin & Shaw (2007). Best fit was obtained for each 
light curve using the same spots and the corresponding parameters are given 
in Table 3. The parameters presented in the last column is the weighted mean 
of the parameters obtained from individual light curves. Coughlin & Shaw 
(2007) analyzed only V- and R-light curves. Our analysis gives very different 
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TABLE 3 

RESULTS OF THE V-, R- AND LBANDPASS LIGHT CURVE ANALYSIS 

FOR NSVS 1186. THE ADOPTED VALUES ARE THE WEIGHTED 

MEANS OF THE VALUES DETERMINED FROM THE INDIVIDUAL 

LIGHT CURVES. 



Parameters 


V 


R 


I 


Adopted 


i° 


87.0±0.5 


88.6±0.2 


88.3±0.2 


88.0±0.4 


Te/A (K) 


5 250[Fix] 


5 250 [Fix] 


5 250 [Fix] 


5 250 [Fix] 


Te//. (K) 


5 040±12 


4 990±14 


5110±34 


5 020±16 


ni 


4.242±0.049 


4.274±0.049 


4.260±0.072 


4.262±0.052 


^2 


3.810±0.026 


3.739±0.020 


3.989±0.054 


3.790±0.026 


ri 


0.289±0.004 


0.288±0.004 


0.286±0.006 


0.288±0.004 


r2 


0.263±0.003 


0.267±0.002 


0.250±0.005 


0.264±0.003 


Li/(Li + L2) 


0.604±0.009 


0.586±0.009 


0.591±0.013 




Y.{w{o^cf) 


0.041 


0.026 


0.131 




weight 


3 


5 


1 





parameters when compared with those obtained by them. They estimate very 
low effective temperatures for the components with respect to their predicted 
masses. The fractional radii for the primary and secondary stars obtained by 
them are 15 % and 20 % smaller than obtained by us using the same light 
curves. In contrary, the absolute radii are 5 % and 3 % larger than we found. 
The main difference is the locations of the spots on the surface of the stars. 
The spot parameters given by them for the stars should be interchanged. One 
spot on the south pole of the primary and one spot on the north pole of the 
secondary well represent the three-color light curves. The spot coverage on 
the surface of the stars is about 8 %. 

3.3. Absolute Parameters 



Combining the parameters obtained by radial velocity and light curve 
analyses we calculated the physical parameters of the components as listed 
in Table 4. We compare the locations of the components in the effective 
temperature- luminosity diagram, i.e. Hertzsprung- Russell diagram (Figure 4). 
While the location of the more massive primary star agrees well with a 
star of having the mass, the secondary component appears to be more lu- 
minous, more probably originating larger effective temperature and radius. 
We also computed the radii of the components as Ri=O.983±O.O3OR0 and 
R2=O.9O1±O.O26R0 for the primary and secondary stars, respectively. 

We computed the luminosities of the components using these radii and 
effective temperatures as Li =0.649±0.069Lo and L2=O.456±O.O49L0 for 
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Fig. 3. The phase folded V- , R- and I-bandpass hght curves for NSVS 1186. The 
best fitting solutions represented by the solid lines are also plotted for comparison 
(see text). 
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TABLE 4 
FUNDAMENTAL PARAMETERS OF THE SYSTEM. 
Parameter Primary Secondary 

Mass (Mo) 0.870±0.074 0.607±0.053 

Radius (Ro) 0.983±0.030 0.901±0.026 

log g (cgs) 4.392±0.019 4.311±0.021 

Teff (K) 5 250±135 5 020±135 

log (L/Lq) -0.188±0.045 -0.341±0.046 

(wsin i)obs. (kms-i) 77±3 62±4 

(i^sin i)caic. (km s^^) 83±3 76±2 



the primary and secondary stars, respectively. We compare the positions of 
the components in L-Te diagram with evolutionary tracks of pre- and post- 
main sequence evolution and as well as with the ZAMS, as shown in Figure 4. 
For constructing solar mctallicity ZAMS and evolutionary tracks, we used the 
Cambridge version of the the STARS code which was originally developed by 
Eggleton (1971) and had substantial revision up to date (Eldridge & Tout, 
2004). The current version is also capable to account molecular contribution 
on equation of state according to Marigo (2002) approximation which may 
be crucial for low mass stars. The less massive component seems to hotter 
and more luminous with respect to the models. The physical parameters 
determined by us are very different from those obtained by Coughlin and 
Shaw (2007). While they estimate effective temperatures between 4 370 and 
5 250 K from the observed colours they adopted an effective temperature of 
3 970 K for the more massive star. The masses of the components estimated 
to be 0.94 and 0.87 and radii of 1.03 and 0.93, and an effective temperature of 
3 750 K for the secondary star. The locations of the component stars in the 
HR diagram correspond to stars having masses below 0.4 M0. The J-,H-, and 
K- band extinctions in the direction of the variable taken from Schlegel et al. 
(1998) as Aj=0.046, Aj:/=0.030 and Ak=0.019 magnitudes. We computed the 
absolute magnitudes as Mj=3.96±0.10, M//=3.50±0.07 and Mx=3.46±0.07 
for the primary star. Then using the JHK magnitudes given in § 2.1.1 we 
estimate the distance to the variable as 668±25, 689±24 and 678±22 parsec 
in J,H, and K bands, respectively. 

4. OVERSIZED STARS IN THE LOW-MASS ECLIPSING BINARIES 

Abt (1963) and Duquennoy, Mayor and Halbwachs (1991) report that bi- 
nary stars are more common than single stars at masses above that of the 
Sun. In contrary, Reid & Gizis (1997) and Delfosse et al. (1999) suggest that 
binaries are not very common in low-mass stars, although approximately 75 
% of all stars in our Galaxy are low-mass dwarfs with masses smaller than 0.7 
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Fig. 4. Comparison between evolutionary models and the physical parameters of 
NSVS 1186 in the effective temperature-luminosity diagram where the luminosi- 
ties determined from the light and radial velocity curves. Theoretical evolutionary 
tracks, both pre- and post- main sequence, for masses of 0.870 M0 (dotted line) 
and 0.607 Mq (dashed line) and ZAMS (solid line) with solar abundance are calcu- 
lated with the STARS code. The circle and square denote the primary and secondary 
stars. 
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Mq. Due to the low binary fractions and their faintness very few low- mass 
eclipsing binary systems have been observed so far both photometrically and 
spectroscopically, yielding accurate physical parameters. Recently in a pio- 
neer study, Ribas (2003,2006) collected the available masses and radii of the 
low-mass stars and compared with those obtained by stellar evolution models. 
Even small numbers of the sample, a total of eight double-lined eclipsing bi- 
naries whose mass and radius are determined with an accuracy of better than 
3 %, the comparison evidently revealed that the observed radii are systemat- 
ically larger than the models. However, the effective temperatures are cooler 
than the theoretical calculations, being the luminosities in agreement with 
those of single stars with the same mass. This discrepancy between the mod- 
els and observations has been explained by Mullan and MacDonald (2001), 
Torres et al. (2006), Ribas (2006) and Lopez- Morales (2007) and others by the 
high level magnetic activity in the low-mass stars. They suggest that stellar 
activity may be responsible for the observed discrepancy through inhibition 
of convection or effects of a significant spot coverage. However, Berger et al. 
(2006) compared the interferometric radii of low-mass stars with those esti- 
mated from theory and find a correlation between an increase in metallicity 
and larger-than-expected radius, i.e. the larger radii are caused by differences 
in metallicity. 

The number of double- lined low-mass eclipsing binaries having precise pho- 
tometric and spectroscopic observations is reached to 21 according to knowl- 
edge of the authors. Combining the results of the analysis of both photometric 
and spectroscopic observations accurate masses, radii, effective temperatures 
and luminosities of the components have been obtained. In Table 5 we list 
absolute parameters for the low-mass stars with their standard deviations. In 
Figure 5 we plot the radii of the low-mass stars in the eclipsing binary systems 
as a function of mass. Theoretical mass-radius (M-R) diagrams for the zero- 
age main-sequence stars with solar abundance taken from the STARS code and 
Siess et al. (2000) are also plotted for comparison. 

The stars below 0.3 M0 appear to trace the theoretical M-R diagram. 
However, larger mass stars, above 0.3 M© , significantly deviate from the the- 
oretical M-R digram. We also computed projected rotational velocities of the 
components using their radii and orbital period. Since the orbital periods are 
generally short we assume spin-orbit synchronization. The rotational veloc- 
ities are plotted in Figure 6 as a function of the radius difference calculated 
as (Robs — Rzams)/ Rzams- This plot shows that as the rotational velocity 
increases, the difference in radius, gets larger for the stars given in Table 5. 

A low-mass star has a convective atmosphere with a radiative core. The 
depth of the convective zone is about 0.28 R for a star of mass about 0.9 
Mq, and gradually increases to about 0.41 R for a star with a mass about 0.4 
Mq. Stars are thought to be fully convective below about 0.35 M©. Therefore, 
the internal structure of such a low-mass star with a deep convective zone is 
tightly dependent on the mixing length parameter, a. In Figure 7a and b we 
compare locations of the low- mass stars in the M-R and M-Tg// diagrams with 
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the theoretical calculations for various a parameters. Convection is modelled 
by mixing length theory (Bohm-Vitense 1958) with the ratio of mixing length 
to pressure scale height (i.e. a = l/Hp). Theoretical models show that there 
is no significant separation in the M-R relation for the stars below about 0.6 
M0 , depending on the mixing length parameter. For masses above this value 
a separation is revealed. However, as the a parameter increases the computed 
effective temperature gets higher. A separation is clearly seen even for the 
masses below 0.6 Mq . We also plotted empirical parameters of the low-mass 
stars in these panels. This comparison clearly exposes that the discrepancies 
in radii and effective temperatures could not be explained only by changing a 
parameter as it is obvious in M-R and M-Tg// plots represented in Figure 7. 
On the other hand Demory et al. (2009) showed that there is no significant 
correlation between metallicity and radius of the single, low-mass stars. An 
alternative explanation remains to be the magnetic activity responsible for 
the observed larger radii but cooler effective temperatures. As demonstrated 
by D'Antona et al. (2000), MuUan and MacDonald (2001) and Chabrier et 
al.(2007) magnetic fields change the evolution of low-mass stars. Due to the 
high magnetic activity in the fast-rotating dwarfs their surface are covered 
by dark spots. Spot coverage in active dwarfs yields larger radii and lower 
effective temperatures. 

5. DISCUSSION 

We obtained the radial velocity curves for the low-mass eclipsing binary 
NSVS 1186. We re-analysed the existing V-, R- and I-bandpass light curves. 
Since the light curves are distorted they are modelled assuming one spot 
located on each star. Then combining the photometric and spectroscopic 
solutions we derived the absolute physical parameters for both components. 
The standard deviations of the parameters have been determined using the 
JKTABSDIM code (Southworth and Clausen 2005). The less massive component 
appears to have a larger radius and higher temperature with respect to the its 
mass. On the other hand the primary star seems to have a larger radius but 
the effective temperature is in agreement with that expected from its mass. 
We suppose that the larger radii of both components may be originated from 
a higher activity in the system at the time that photometry is performed. 
The locations of the components in the H-R diagram have been plotted and 
compared with theoretical models. The secondary component seems to have 
slightly brighter with respect to its mass. 

Basic physical parameters determined from double lined mass binaries , 
such as mass, radius, effective temperature and orbital period, for low-mass 
stars are collected. The locations of the low-mass stars in the mass-radius 
panel are compared with the theoretical ZAMS models. Since the low-mass 
stars evolve very slowly, we prefer a comparison with the ZAMS models. The- 
oretical models tend to underestimate the radii at a given mass in the range 
of 0.3 and 0.6 Mq . Due to the evolutionary effect on the radius, we take 
into account only this mass range, preventing any confusion. Since low-mass 
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TABLE 5 

THE MASS, RADIUS, EFFECTIVE TEMPERATURE AND ORBITAL 

PERIOD OF THE LOW-MASS STARS IN THE DOUBLE-LINED 

ECLIPSING BINARIES. 



Object 



Star 



Mass/M0 



Radius/RQ 



Tcff/K Porb/day Refs 



CM Dra 


A 




B 


YY Gem 


A 




B 


CU Cnc 


A 




B 


GU Boo 


A 




B 


TrES-HcrO-07621 


A 




B 


2MASS J05162881-I-2607387 


A 




B 


UNSW-TR 2 


A 




B 


NSVS 06507557 


A 




B 


NSVS 02502726 


A 




B 


T-Lyrl-17236 


A 




B 


2MASS J01542930+0053266 


A 




B 


GJ 3236 


A 




B 


SDSS-MEB-1 


A 




B 


BD -22 5866 


A 




B 


NSVS 01031772 


A 




B 


NSVS 11868841 


A 




B 


GJ 2069A 


A 




B 


2MASS J04463285+1901432 


A 




B 


NSVS 6550671 


A 




B 


IM Vir 


A 




B 


RXJ0239.1 


A 



0.2310 ± 0.0009 

0.2141 ± 0.0008 

0.5992 ± 0.0047 

0.5992 ± 0.0047 

0.4349 ± 0.0012 

0.3992 ± 0.0009 

0.6101 ± 0.0064 

0.5995 ± 0.0064 

0.4930 ± 0.0030 

0.4890 ± 0.0030 

0.7870 ± 0.0120 

0.7700 ± 0.0090 

0.5290 ± 0.0350 

0.5120 ± 0.0350 

0.6560 ± 0.0860 

0.2790 ± 0.0450 

0.7140 ± 0.0190 

0.3470 ± 0.0120 

0.6795 ± 0.0107 

0.5226 ± 0.0061 

0.6590 ± 0.0310 

0.6190 ± 0.0280 

0.3760 ± 0.0170 

0.2810 ± 0.0150 

0.2720 ± 0.0200 

0.2400 ± 0.0220 

0.5881 ± 0.0029 

0.5881 ± 0.0029 

0.5428 ± 0.0027 

0.4982 ± 0.0025 

0.8700 ± 0.0740 

0.6070 ± 0.0530 

0.4329 ± 0.0018 

0.3975 ± 0.0015 

0.4700 ± 0.0500 

0.1900 ± 0.0200 

0.5100 ± 0.0200 

0.2600 ± 0.0200 

0.9810 ± 0.0120 

0.6644 ± 0.0048 

0.7300 ± 0.0090 

0.6930 ± 0.0060 



0.2534 ± 0.0019 

0.2398 ± 0.0018 

0.6194 ± 0.0057 

0.6194 ± 0.0057 

0.4323 ± 0.0055 

0.3916 ± 0.0094 

0.6270 ± 0.0160 

0.6240 ± 0.0160 

0.4530 ± 0.0600 

0.4520 ± 0.0500 

0.7880 ± 0.0150 

0.8170 ± 0.0100 

0.6410 ± 0.0500 

0.6080 ± 0.0600 

0.6000 ± 0.0300 

0.4420 ± 0.0240 

0.6740 ± 0.0600 

0.7630 ± 0.0050 

0.6340 ± 0.0430 

0.5250 ± 0.0520 

0.6390 ± 0.0830 

0.6100 ± 0.0930 

0.3828 ± 0.0072 

0.2992 ± 0.0075 

0.2680 ± 0.0090 

0.2480 ± 0.0080 

0.6140 ± 0.0450 

0.5980 ± 0.0450 

0.5260 ± 0.0028 

0.5088 ± 0.0030 

0.9830 ± 0.0300 

0.9010 ± 0.0260 

0.4900 ± 0.0800 

0.3300 ± 0.0400 

0.5700 ± 0.0200 

0.2100 ± 0.0100 

0.5500 ± 0.0100 

0.2900 ± 0.0100 

1.0610 ± 0.0160 

0.6810 ± 0.0130 

0.7410 ± 0.0040 

0.7030 ± 0.0020 



3130 


± 


70 


1.276 


1 


3120 


± 


70 






3820 


± 


100 


0.820 


2 


3820 


± 


100 






3160 


± 


150 


2.794 


3 


3125 


± 


150 






3920 


± 


130 


0.492 


4 


3810 


± 


130 






3500 


± 




1.137 


5 


3395 


± 








4200 


± 




2.619 


6 


4154 


± 










± 




2.144 


7 




± 








3960 


± 


80 


0.520 


8 


3365 


± 


80 






4300 


± 


200 


0.560 


9 


3620 


± 


205 






4150 


± 




8.430 


10 


3700 


± 








3730 


± 


100 


2.639 


11 


3532 


± 


100 






3310 


± 


no 


0.770 


12 


3240 


± 


no 






3320 


± 


130 


0.410 


13 


3300 


± 


130 








± 




2.211 


14 




± 








3615 


± 


72 


0.368 


15 


3513 


± 


31 






5260 


± 


no 


0.602 


This paper 


5020 


± 


no 








± 




2.771 


16 




± 








3320 


± 


150 


0.630 


18 


2910 


± 


150 






3730 


± 


60 


0.193 


19 


3120 


± 


65 






5570 


± 


100 


1.309 


20 


4250 


± 


130 






4645 


± 


20 


2.072 


21 


4275 


± 


15 







Ref: (l)Morales et al. (2009a), (2)Torres & Ribas (2002), (3)Ribas (2003), (4)L6poz- 
Morales & Ribas (2005), (5)Creevey et al. (2005), (6)Bayless & Orosz (2006), (7)Young 
et al. (2006), (8)Qakirli & Ibanoglu (2010), (9)gakirli, Ibanoglu & Giingor (2009), (10) 
Devor et al. (2008), (ll)Becker et al. (2008), (12)Irwin et al. (2009), (13)Blake et al. 
(2008), (14)Shkolnik et al. (2008), (15)L6pez-Morales et al. (2006), (16)Delfosse et al. 
(1999), (17)Beatty et al. (2007), (18)Hebb et al. (2006), (19)Dimitrov & Kjurkchieva 
(2010) ,(20)Morales et al. (2009b), (21)L6pez-Morales fc Shaw (2007) 
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K 0.6 - 



0.4 0.6 



Fig. 5. Locations of the low-mass stars in the mass-radius diagram. The components 
of NSVS 1186 are shown by triangles. The solid and dashed lines represent the ZAMS 
models taken from the STARS code and Siess et al.(2000), respectively, with solar 
abundance. 
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Fig. 6. The rotational velocities are plotted against the radius difference from the 
ZAMS with the same mass. 



stars have deeper convective layers in their atmospheres we also computed 
the radii and temperatures for the above-mentioned mass range for various 
mixing-length parameters. We find that mixing-length parameter does not 
significantly affect the radii, in contrary, it changes considerably the effective 
temperatures of the stars. Comparison with the observations shows that the 
observed larger radii of the late K and M dwarfs can not only be represented 
with the values of mixing-length parameters. These comparisons confirm that 
magnetic field induced inhibition of convection in fast-rotating low-mass stars 
with deep convective zones lead to larger radii and lower effective temper- 
atures than those expected from theoretical models (D'Antona et al. 2000, 
MuUan and MacDonald 2001, Chabrier et al. 2007) . 
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Fig. 7. Mass-radius (a) and mass-effective temperature (b) relationships for zero- 
age main-sequence low-mass stars calculated with the STARS code for values of the 
mixing length parameter as q = 4 (dotted line), a = 2 (solid line), a = 1 (dashed 
line) and q = 0.5 (dot-dashed line). The values of the alpha parameter higher 
than 0.5 appear to effect the radius slightly for a star below 0.7 Mq. As the alpha 
increases the effective temperatures get higher. 
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